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O ' Abstract 

\q . Following a suggestion given in [1], we show how the U(1)®Z 2 symmetry of 

the fully frustrated XY (FFXY) model on a square lattice can be accounted for 
in the framework of the m-reduction procedure developed for a Quantum Hall 
system at "paired states" fillings v — 1 [2, 3]. The resulting twisted conformal 
field theory (CFT) with central charge c = 2 is shown to well describe the 
Ph. physical properties of the FFXY model. In particular the whole phase diagram 

is recovered by analyzing the flow from the Zi degenerate vacuum of the c = 2 
CFT to the infrared fixed point unique vacuum of the c = | CFT. The last 
theory is known to successfully describe the critical behavior of the system at 
the overlap temperature for the Ising and vortex-unbinding transitions. 
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1 Introduction 



A few years ago Foda [1] proposed that the continuum limit of a FFXY model is equivalent, 
at the overlap temperature, to a free massless field theory of one boson and one Majorana 
fermion, that is a superconformal field theory with central charge c = 3/2. Such a model 
shares both discrete and continuous symmetries, which describe well a Josephson junction 
array in a transverse magnetic field, provided that the external magnetic flux threading 
each cell of the array is |<E»o, where $ — ff is the superconducting flux quantum [4, 5]. 

The XY model on a square lattice in the presence of an external magnetic field transver- 
sal to the lattice plane is described by the action: 

H = ^ C ° S ^ ~ ~ Aij ^ ' ^ 

{ij) 

where {<p} are the phase variables defined on the sites, the sum is over nearest neighbors, 
J > is the coupling constant and A; L j — ^ J? A-dl is the line integral along the bond 
between adjacent sites i and j. We consider the case where the bond variables Ay are fixed, 
uniformly quenched, out of equilibrium with the site variables and satisfy the condition 
^2 p Aij = 27r/; here the sum is over each set of bonds of an elementary plaquette and 
/ is the strength of frustration. We assume that the local magnetic field in Eq. (1) is 
equal to the uniform applied field; such an approximation is more valid the smaller is the 
sample size L compared with the transverse penetration depth Aj_ [6]. In the case of full 
frustration, i.e. / = |, of interest to us here, such a model has a continuous U (1) symmetry 
associated with the rotation of spins and an extra discrete symmetry, as it has been 
shown analysing the degeneracy of the ground state [7, 8]. Choosing the Landau gauge, 
such that the vector potential vanishes on all horizontal bonds and on alternating vertical 
bonds, we get a lattice where each plaquette displays one antiferromagnetic and three 
ferromagnetic bonds. Such a choice corresponds to switching the sign of the interaction 
term in Eq. (1) and is closely related to the presence of two ground states with opposite 
chiralities, the first one invariant under shifts by two lattice spacings and the second one 
invariant under shifts by one lattice spacing. 

In the presence of such a degeneracy we can use the m-reduction technique [3] which 
has been successfully applied to a quantum Hall fluid in [2, 3, 9, 10, 11] and to a fully 
frustrated Josephson junction ladder with Mobius boundary conditions in [12]. In Section 
4 we will see how the discrete version of the m-reduction procedure, successfully used to 
build a CFT with a Z m -twist, well describes such a situation. 

The Zi symmetry of the FFXY model is broken at low temperature and will be restored 
beyond a certain temperature after the formation of domain walls separating islands of 
opposite chirality. The Ising transition overlaps to a vortex- unbinding transition [13], 
which is associated with the U(l) symmetry [8, 14]. 

The action (1) can be cast into a form where both the U(l) and Z2 symmetries are 
manifest, through the Villain approximation [15]. In this way the spin- wave and the 
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vortex contributions can be separated. Furthermore, by integrating out the spin waves, 
the resulting vortex contribution can be rewritten as a fractionally charged Coulomb gas 
(CG) defined on the dual lattice [8, 16]: 

H = -^fJ2( m ( r ) + /) G ^ r ') ( m ( r ') + /) • ( 2 ) 

Here we have lim| r _ r /|_ +00 G (r, r') = log \r — r'\+^n and the neutrality condition ^ r (m (r) + /) = 
Yl r n(r) = must be satisfied. It is now evident that the ground state for / = | consists 
of an alternating lattice of logarithmically interacting ±| charges and is doubly degener- 
ate. Such a model exhibits two possible phase transitions, an Ising and a vortex-unbinding 
one, and their relative order has been deeply studied in the literature [7]. 

It is just the doubly degenerate checkerboard pattern of vortices in the ground state 
that gives rise to the Ising-like Z 2 discrete symmetry in addition to the Kosterlitz-Thouless 
(KT) continuous U(l) one, associated with the uniform rotation of all phase angles [13]. 
The issue whether there are two distinct phase transitions, Ty > Tjsing or Ty < Tjsing 
with T v and T ISING marking respectively the breaking of U(l) and of Z 2 symmetry 
[17], or a single transition with the simultaneous breaking of both symmetries [18] has 
been widely investigated. Indeed an exact mapping of the classical 2D FFXY model 
onto an alternating 19- vertex model has been established [19], whose equivalence with the 
restricted solid-on-solid (RSOS) model is well known [20]. In turn the RSOS model can be 
viewed as a discretized time generalization of the spin-1 one dimensional quantum chain, as 
shown in detail in Ref.[21] by means of the transfer matrix formalism. The phase diagram 
of such a system has been widely investigated in a simplified case, i.e. the XXZ quantum 
chain with single-ion-type anisotropy, by means of analytical and numerical techniques 
[21, 22, 23]. It shows both U(l) and Z 2 symmetry breaking in the space D — J z (where the 
parameter D is the uniaxial single-ion anisotropy) [21, 24, 25]. Depending on the values 
of the parameters different phases of the system emerge: the Haldane phase, the large-/} 
phase, two XY phases (XY1 and XY2) with different Z 2 symmetry, the ferromagnetic 
phase and the Neel phase. Then various types of phase transitions take place. In particular 
all the three possibilities above quoted can be recognized: a gapful-gapless KT transition 
between the Haldane phase and the XY\ one followed by a transition between the two 
different XY phases, XY1 and XY2; an Ising transition between the Haldane phase and 
Neel phase and after that a transition Neel-XF2 at J z = for large negative D; an XY 
and Ising transition lines which merge at a tricritical point within c = | CFT universality 
class. All that reveals the richness of the phase structure of a class of frustrated systems 
with the same ground state but the controversy on the phase diagram structure of the 
FFXY model remains still unsolved. 

Recently Korshunov [26] argued that a new transition could take place in the 2D FFXY 
model, well below the bulk transition, due to the unbinding of kink-antikink pairs on the 
domain walls associated with Z 2 symmetry. Such a transition could lead to a decoupling 
of phase coherence across domain boundaries, so producing two distinct bulk transitions 
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with T v < Tjsing and then a coupled XY-Ising model [27, 28]. Direct numerical evidence 
supporting the existence of such a transition at a temperature T w below the bulk ones 
has by now been provided [29]. Also Teitel and Jayaprakash [14] have considered the 
case Ty < Tising in terms of a dual Coulomb gas model. They argued that the helicity 
modulus Y would decrease from its T = value by means of fluctuations producing dipole 
moments. So, if Ising domains of size Rising carried a total dipole moment proportional 
to their size, they would make Y = in a continuous way, due to the divergence of 
Rising as T — > Tjsing- At T = it has been found though that the only Ising-like 
domains are those which carry no dipole moment: they cannot produce a reduction in 
Y. But at T > Tw, according to Ref. [26], kink-antikink excitations on the boundaries 
of the domain walls unbind, so the Ising domains would acquire large dipole moments 
without cost in free energy. Such a conclusion [29] supports strongly the existence of two 
separate transitions with Ty < Tising- This scenario has been recently confirmed also by 
numerical simulations on very large lattices and by means of a careful finite-size scaling 
analysis in Refs. [30] . The critical behaviour at the overlap temperature is well described 
by a conformal field theory with central charge c = 3/2 [1]. The interesting scenario 
just outlined is the result of an interplay between the Ising and gaussian sectors of the 
model: they are linked through the fractional vortices which reside at the corners of the 
domain walls in the off-critical theory at the microscopic level [7]. More recent Monte 
Carlo simulations of 2D XY-type models [31] appear to confirm the above scenario [32]. 
Indeed the sequence of phase transitions on going from higher to lower temperatures is: 
first Ising, then KT and next kink-antikink unbinding. This picture will then remain 
(with Ising and KT transition very close to each other) until a particular value of the 
parameters where all the three lines merge at the same time. After that the transition is 
first order and the critical point is characterized by a central charge c = 2 [32]. 

The FFXY model has an interesting experimental realization as a Josephson junctions 
array (JJA) in a transverse magnetic field with half a magnetic flux quantum per cell 
[33, 34] . Measurements have been performed on the electrical properties of such a system, 
in particular on the resistance as a function of the temperature and the external magnetic 
field which is closely related to the helicity modulus and then to the vortex-unbinding 
transition [33]. The fact that the vortex transition is shifted downwards in temperature 
[33] is in agreement with the results found in Ref. [35] , which refer to an XY model with 
a modified potential of the form: 

V (</?) = cosy? + acos2(/9. (3) 

That tentatively leads to the conclusion that the arrays of Ref. [33] display all the 
properties of the FFXY model with higher harmonics contributions to the potential [1]. 

In this paper we construct a CFT for such systems which extends the results of Ref. 
[1], so recovering the whole phase diagram quoted in Refs. [27, 28, 36, 37], then we 
close with a brief account on the relation between the FFXY model and the physics of a 
Josephson junction array in the classical approximation [38]. 
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The paper is organized as follows. 

In Section 2 we discuss in a qualitative way the phase diagram of the FFXY model 
describing some models which have the same U(l) <g) Z 2 degenerate ground state and are 
believed to be in the same universality class. 

In Section 3 we recall some aspects of the m-reduction procedure, in particular we 
show how the m = 2, p = case well accounts for the U(l) (g> Z 2 symmetry of the FFXY 
model. In such a framework we give the whole primary fields content of the theory on 
the plane. 

In Section 4 we give a discrete version of the m-reduction procedure and then, starting 
from our CFT results, we show how to recover the whole FFXY phase diagram introduced 
in Section 2. 

In Section 5 some comments and outlooks are given. 



2 Phase diagram of the FFXY model 



In the literature the phase transitions of the FFXY model have been studied analyzing 
other models which have the same U(l) <E> Z 2 degenerate ground state and are believed 
to be in the same universality class. In particular, by using symmetry arguments[7, 18] 
or a Hubbard-Stratonovich transformation [27], the FFXY model can be reformulated in 
terms of a system of two coupled XY models with a symmetry breaking term: 



H = A 



E E cos (^0-^(0) 

i=i,2 ( r y) 



+ hJ2cos2(^ 1 \r)- V (2 \r)) . (4) 



Such a system has been studied in detail by Granato and Kosterlitz [39]. The limit 
h — > corresponds to a full decoupling of the fields ip( l \ i — 1, 2, so giving rise to a CFT 
with central charge c = 2 which describes two independent classical XY models (in the 
continuum limit). In the h — > oo limit, the two phases ip^\ i — 1, 2 are locked [39]: 



(p W(r)-<pW(r)=irj, j = 1,2; 



,(2) 



(5) 



as a consequence the model gains a symmetry U (1) ®Z 2 and its Hamiltonian renormalizes 
towards a model described by: 



H = H (h -> oo) = A (! + s rV) cos (V 1} (0 - 

( r i r ' ) 



(6) 
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where (p^(r) and s r = cosnj = ±1 are planar and Ising spins respectively. In this way a 
model is obtained which is consistent with the required symmetry, the XY-Ising one and 
whose Hamiltonian has the general form [27, 36] : 



(1 + s r v) cos (y? (1) (r) - y? ( V)) + Cs r s r , 



(7) 



( 




which is believed to be in the universality class of the classical FFXY model. However 
such a model does not contain fractional vortices at the corners of the domain walls [7], 
which are believed to give rise to an interesting interplay between the continuous and 
discrete symmetries [1]. It is just this complex interplay which accounts for the unusual 
and not yet fully understood critical behaviour of the FFXY model on the square lattice 
[1, 27, 28, 36, 37], that will be discussed in the following. 

The phase diagram of the model [27, 28, 36, 37] with Hamiltonian (7), as depicted in 
Fig. 1, is built up with three branches which meet at a multi-critical point P. The branch 
PT corresponds to single transitions with simultaneous loss of XY and Ising order while 
the other two describe separate Kosterlitz-Thouless and Ising transitions. The critical 
line PT becomes a first order one at the tricritical point P; it seems to be non-universal 
[27, 40] as the values of the critical exponents {y = 0.79, rj = 0.40) estimated by finite- 
size scaling of large systems are inconsistent with pure Ising critical behavior [y = 1, 
V — \)- Also the numerical estimate for the central charge c (c ~ 1.60) [28] is higher 
than the expected value (c = 3/2) for a critical behavior given by the superposition of 
critical Ising (c = 1/2) and gaussian (c = 1) models [1]. Indeed the central charge seems 
to vary continuously from c 1.5 near P to c 2 at T [28]. All such results would 
suggest the existence of a parameter changing along the critical PT line that does not 
affect the symmetry [37]. The system seems to be not conformally invariant [27], so it 
would be natural to consider all possible integrable perturbations of a superposition of 
critical Ising and XY models as a starting point to study the vicinity of the point P [37]. 
Furthermore different numerical techniques have led to different results suggesting either 
single [18] or double transitions [17] along the critical PT line. How to reconcile such 
conflicting results? A possible solution is suggested in Ref. [41] where the Coulomb gas 
Hamiltonian of Eq. (2) gets modified through an additional antiferromagnetic coupling 
between nearest-neighbour vortices: 



Such a model has a phase diagram with a structure close to the XY-Ising model one 
of Fig. 1 with the occurrence of double or single transitions on the PT line depending on 
this additional coupling; for J ^ the two transitions separate with the vortex-unbinding 
one occurring at a lower temperature [41]. A similar effect has been observed theoretically 
in Ref. [35], where a XY model with various modifications of the cosine- type potential 
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is studied, and experimentally in Ref. [33]: it can be ascribed to the presence of higher 
harmonics contributions to the potential. Such a finding is a crucial one, as we will show 
in Section 4. 

Another route to the solution of all these issues is the m-reduction technique [2, 3] 
which describes well the two fixed points with c = 3/2 and c = 2 and accounts for 
non trivial boundary conditions for the square lattice. Such a realization of the FFXY 
model on closed geometries could be relevant for the description of JJAs with non trivial 
topologies, which are believed to provide a physical implementation of an ideal quan- 
tum computer [42] because of the topological ground state degeneracy "protected" from 
external perturbations [43]. We introduce the m-reduction in the next Section. 

3 A review of the m-reduction procedure 

In this Section we recall those aspects of the twisted model (TM) which are relevant for 
the FFXY model. We focus in particular on the m-reduction procedure for the special 
m = 2 case (see Ref. [2] for the general case), since we are interested in a system with 
U(l) ®Z 2 symmetry. We showed in Refs. [2, 3] that such a theory describes well a system 
consisting of two parallel layers of 2D electron gas in a strong perpendicular magnetic 
field, with filling factor = f° r eacn °f the two a = 1, 2 layers and total filling 
v = z/ 1 ) + z/ 2 ) = ^j-j-. Regarding the integer p, characterizing the flux attached to the 
particles, we choose the "bosonic" value p — 0, since it enables us to describe the highly 
correlated system of vortices with flux quanta 

Let us start from the bosonic "filling" v = |, described by a CFT with c = 1 in terms 
of a scalar chiral field Q compactified on a circle with radius R 2 = l/u = 2. It is explicitly 
given by: 



with a n , q and p satisfying the commutation relations [a n , a n i] = n5 n ,n' and [q,p] = %■ 

From such a CFT (mother theory), using the m-reduction procedure, which consists in 
considering the subalgebra generated only by the modes in Eq. (9) which are a multiple 
of an integer m, we get a c = m orbifold CFT (daughter theory, i.e. the TM) which, 
for m = 2, in the next Section will be shown to describe the whole phase diagram of 
the FFXY model. Then the fields in the mother CFT can be organized into components 
which have well defined transformation properties under the discrete Z m (twist) group, 
which is a symmetry of the TM. By using the mapping z — > z 1 !'" 1 and by making the 
identifications a nm+ i — ► y/ma n+ i/ m , q — > -7=?, the c = m CFT (daughter theory) is 




(9) 



obtained. 
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Its primary fields content, for the special m = 2 case, can be expressed in terms of a 
Z 2 -invariant scalar field X(z), given by 

X(z) = ±(Q(z) + Q(-z)), (10) 
describing the electrically "charged" sector of the new theory, and a twisted field 

*{z) = \lQ{z)-Q{-z)), (11) 

which satisfies the twisted boundary conditions 0(e M z) = —(/>(z) and describes the "neu- 
tral" sector [2]. 

The TM primary fields are composite vertex operators V (z) = U x (z) ip (z) where 

U x (*) = "7= : e iaX{z) : , o? = 2 (12) 

V z 

is the vertex of the "charged" sector and 

V> (*) = -7= : e ta4>(z) ■ (13) 

V z 

is the "neutral" one. We point out that in the neutral sector it is useful to introduce the 
two chiral operators: 

i) (z) = (: e ia *W : + : e <"*(-*) .) ; ^(z) = (: e^ (z) : - : e ia ^ z) :) ; 

the first one does not change the boundary conditions while the second one does. In a 
fermionized version of the theory they correspond to two c = 1/2 Majorana fermions with 
Ramond and Neveu-Schwartz boundary conditions [2, 3]. Furthermore in the TM they 
appear to be not completely equivalent. 

In fact the whole TM theory decomposes into a tensor product of two CFTs, a twisted 
invariant one with c = 3/2 (the Moore-Read (MR) theory with symmetry U(l) ® Z 2 ) and 
the remaining c = 1/2 one realized by a Majorana fermion in the twisted sector. Such a 
factorization can be unambiguously pointed out on the torus topology where the partition 
function of the c = 2 theory, i.e. the TM, can be written as 5 : 

Z(w c \t) = Z mr (w c \t)Z 7 (t). (15) 

Here Z MR (w c \r) stands for the MR partition function and Zj (r) is the partition function 
describing the Ising degrees of freedom which are Z 2 antisymmetric (i.e. ip (z) of Eq. (14) 
in the plane geometry). 

5 See Appendix for the explicit expression and Ref. [3] for its derivation. 
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Let us point out that the energy-momentum tensor of the Ramond part of the neutral 
sector develops a cosine term: 



a clear signature of a tunneling phenomenon which selects out a new stable vacuum, the 
c = 3/2 one. We identify such a Z 2 invariant theory with the one describing the FFXY 
model conjectured in Ref. [1]. Here we only point out that the critical behavior of the 
c = 3/2 theory will be a superposition of the critical behavior of its components with c = 1 
and c = 1/2 respectively The gaussian component will lead to a gaussian critical line 
that ends at a vortex-unbinding transition point, while the Ising component will lead to 
an Ising critical point. The conformal field theory, considered in this paper, is not meant 
to predict the position of the Ising point with respect to the gaussian critical line, it only 
tells us that there exists an Ising transition overlapping somewhere with a gaussian critical 
line. In the following we exhibit all the relevant physical consequences of our theory. 

4 TM description of the FFXY phase diagram 

In this Section we will derive the phase diagram of the FFXY model 6 in terms of the RG 
flow which originates from perturbing our TM model. 

As we pointed out in Section 2, the limit h — > in the Hamiltonian of Eq. (4) 
corresponds to a full decoupling of the fields (p( % \ i = 1,2. In the continuum limit that 
gives rise to a CFT with two scalar Fubini boson fields ip® and with central charge c = 2. 
Let us stress that a good candidate to describe the FFXY model at criticality around 
the point T of the phase diagram is a CFT, with c = 2, which accounts for the full 
spectrum of excitations of the model. So it is worth recalling briefly the structure of such 
excitations for the FFXY model. A domain wall is a topological excitation of the double 
degenerate ground state and it can be defined as a line of links, each one separating two 
plaquettes with the same chirality. So, through a domain wall the alternating structure 
(the checkerboard pattern) of the ground state is lost. Kinks and antikinks are excitations 
which live on the domain walls and are described by fractional vortices with +1/2 and 
— 1/2 topological charge. Such particle-like excitations belong also to the spectrum of 
the fully frustrated ladders of Josephson junctions and in [12] we have shown how the 
TM model exactly describes them. In particular these excitations can be generated from 
the ground state by closing the ladder, by imposing the coincidence of opposite sides. In 
fact, in this closed geometry, for an even number of plaquettes we obtain the ground state 
while for the odd case this point-like excitation is generated on the boundary. In the TM 
model they are described by different boundary conditions, untwisted and twisted ones, 
respectively. 

6 As it has been described in Section 2 and inferred by Migdal-Kadanoff renormalization [36], Monte 
Carlo simulations [27] and Monte Carlo transfer matrix calculations [28]. 




(16) 
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Similar considerations apply to the FFXY array; in particular the central role played 
by the closed geometry and by the boundary conditions in the description of the excitation 
spectrum is clear already at the level of the lattice model. By imposing the coincidence 
between opposite sides of the square lattice, we obtain a closed geometry, which is the 
discretized analogue of a torus. Now for the ground state two topologically inequivalent 
circumstances arise, one for even and the other one for odd number of plaquettes. In the 
even case the end plaquettes on the opposite sides of the lattice have opposite chirality, 
while in the odd case they have the same chirality. So the ground state on the square 
lattice maps into the ground state for the even case while it generates two straight domain 
walls along the two cycles of the torus for the odd case. Such a behaviour has to be taken 
into account by opportune boundary conditions on the field (p® at the borders of the finite 
lattice. These non trivial boundary conditions naturally arise when we implement the m- 
reduction procedure in the discrete case. To this aim, let (— L/2,0), (L/2,0), (L/2,L), 
(— L/2,L) be the corners of the square lattice C and assume that the fields </?W satisfy 
the following boundary conditions: 

(pW(r) =</? (2) (r) forre£nx, (17) 

where x is the x axis. The above boundary conditions allow us to consider the two fields 
ip^ and (p^ on the square lattice C as the folding of a single field Q, defined on the 
lattice Cq with corners (— L/2, — L), (L/2, —L), (L/2,L), (— L/2,L). More precisely we 
define the field Q as: 

^ [r) ~ \ y?( 2 )(-r) for r G (-£) fl C . ^ 

We can implement now a discrete version of the m-reduction procedure (m = 2) by 
defining the fields: 

X(r) = ±(Q(r) + Q(-r)), (19) 



S(r) = -(Q(r)-e(-r)), 



(20) 



where r E Cq. The resemblance with the continuum version of the m-reduction procedure 
is evident and the fields X and $ are symmetric and antisymmetric with respect to the 
action of the generator g : r —> — r of the discrete group Z 2 . The Hamiltonian in Eq. (4) 
can be rewritten in terms of these fields and, for h — 0, it becomes: 



H = 2A cos ( X ( r ) ~ X ( r ')) cos ( $ ( r ) - $ ( r ')) > 

(r/)ec 

which in the continuum limit corresponds to the action of our TM model: 



(21) 



A 



\ {OXf + \ (9$) 2 



d x. 



(22) 
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It is worth pointing out that the fields X and $ are scalar fields and so the chiral fields 
defined by Eqs. (10), (11) can be seen as their chiral components. Moreover the group 
Z 2 is a discrete symmetry group, indeed both H and A are invariant under its action. 

In the following we will show how the continuum version of the FFXY model and its 
phase diagram can be described by the action: 

d 2 x, (23) 

defined in terms of m = 2 reduced fields X and $, which embodies in its "neutral" 
sector the higher harmonic potential term conjectured in Refs. [1, 35, 41]. We assume 
the constraints (3 2 < 87r, which characterizes both the cosine terms as relevant scalar 
fields perturbing the gaussian fixed point, and \8\ < n/2 (see [44]). Thus the "neutral" 
sector is a two-frequency sine-Gordon theory that can be viewed as a deformation of a 
pure sine-Gordon one with the perturbing term A cos {(3Q/2 + 5). The ultraviolet (UV) 
fixed point fi = 0, A = of the action (23) corresponds to the TM model with central 
charge c = 2, describing the fixed point T in the phase diagram of the FFXY model. 
Here the argument described in [44] can be adapted to our particular case to study the 
RG flow in the "neutral" sector of our theory. Let us define the dimensionless variable 
rj = A /i -( 87r -( /3 / 2 ) 2 )/( 8w - /32 ) ; when rj = (A = 0) the "neutral" sector reduces to a sine- 
Gordon model with a particle spectrum which consists of solitons and antisolitons and, for 
j3 2 < 47r, some breathers. Switching on the perturbation (i.e. for A ^ 0) a confinement of 
solitons into states with zero topological charge takes place and packets formed by 2 of the 
original solitons or antisolitons survive as stable excitations for generic values of | <5 1 < n/2. 
In the limit rj — > oo the 2-soliton evolves into the 1-soliton of the pure sine-Gordon model 
with fi = 0. An unbinding phenomenon takes place in the particular 5 = ±7r/2 case for 
finite f] and the 2-soliton decomposes into a sequence of two kinks K\. In such a case the 
limit rj — > oo implies a transmutation of a composite topological excitation (the two kinks 
K\) into an elementary one (the 1-soliton). So the existence of an intermediate critical 
value f] = r] c is required at which a phase transition takes place and the RG flow ends into 
the infrared (IR) fixed point described by a CFT with central charge c = 1/2, the Ising 
model. The central charge of the full model (23) so changes from c = 2 of the UV fixed 
point to c = 3/2 of the IR fixed point, i.e. we recover early known results of Monte Carlo 
simulations [28]. Indeed in the following we clarify how such an IR fixed point coincides 
with the U(1)®Z 2 symmetric component (c = 3/2) of the whole TM model, which results 
then to properly describe the FFXY model conjectured in Ref. [1], i.e. the fixed point 
P in the phase diagram of the FFXY model. Then the full phase diagram of the FFXY 
model in Fig. 1 would be recovered. 

To such an extent we define fermion fields on two parallel layers in terms of the 
"neutral" boson fields. In particular let ipi (x) be the fermion on layer i, whose chiral 
components at the UV fixed point are defined by: 

m L/R = ± : : ; (24) 
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A 



\ {dX? 



\m 2 



H cos (/3<E>) + A cos I ^$ + 5 



where (&) L / R are the left and the right components of <E>. In terms of them the Z 2 
symmetric ip and Z 2 antisymmetric ip fermions read as: 

^ = ^1 + ^2) , ^=^1-^2). (25) 

The action (23) at \5\ = n/2 can then be rephrased in terms of two interacting Ising layers 
as: 

A = J X - {dX) 2 d 2 x + A[ sing + A{ sina + p J d 2 x(e 1 (x) e 2 (x)) + (26) 

+fi J d 2 x(e 1 (x) + e 2 (x)) + A J d 2 x(a 1 (x)a 2 (x)), (27) 

where A\ s%n9 denotes the fixed point action of the Majorana fermion ipi and £j (x) , 0, (x) 
are the energy and spin of the i-th Ising model (i = 1,2). This correspondence comes 
out observing that the two critical Ising models plus the interaction term in the coupling 
p define a line of c = 1 CFTs parameterized by p since the operator s\s 2 is marginal 
(A £l£2 = 1). In this one-parameter family of c = 1 CFTs the thermal e± + e 2 and the 
magnetic 0\o 2 operators are known [45] to have p-dependent conformal dimensions which 
however satisfy the relation: 

A aia2 (p) A CT1CT2 (0) 1 



(P) (0) 4" 



(28) 



It allows us to say that the thermal and magnetic terms in (27) play the role of the two 
cosine terms in (23), in particular the equivalence works under the following identifications: 

£1 + e 2 ~ cos (/?$) , C7i(72 ~ sin (t^^J ' ( 29 ) 

where the parameters p and (3 are linked by the relation A £l+£2 (p) = f3 2 /87r. Finally we 
are in the position to identify unambiguously the IR fixed point c = 3/2 of our theory. 
Indeed for each value of p(/3) there is a critical line in the /iA-plane corresponding to the 
RG flow of our theory from the UV fixed point c = 2 to the IR fixed point c = 3/2; in 
particular we can consider the RG flow selected by the limit A — > 00. In such a case the 
magnetic term in (27) forces the spins o\ and o 2 to align everywhere reducing the two 
Ising system to a single Ising one. Thus at the IR fixed point the degrees of freedom of the 
fermion ip decouple and we are left with the U(l) <E> Z 2 symmetric component (c = 3/2) 
of the whole TM model, i.e. the MR theory with partition function Z MR (u> c |r), as it can 
be immediately seen from Eq. (15). 



5 Conclusions and outlooks 

In this paper we proposed a c = 2 CFT description of the FFXY model, which extends 
an early proposal by Foda [1]. A crucial role in such a construction is played by the Z 2 
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(Z m in general) discrete symmetry built-in in the 2 (m in general)-reduction technique, 
which allows for non trivial topological properties of the vacuum. 

By perturbing the c = 2 CFT with relevant operators, which provide higher harmonics 
contribution to the potential, it was shown how the unbinding of kink-antikink states can 
give rise to a massless line flow, which fully reproduces the phase diagram for the FFXY 
model. As a result of the flow the non invariant Z 2 degrees of freedom decouple, the 
partition function gets reduced to the Moore-Read one and the degeneracy of the ground 
state lost. The attractive infrared fixed point is identified with the c = § superconformal 
field theory, that is a free massless field theory of one boson and one Majorana fermion, 
which well describes the critical behavior of the FFXY at the overlap temperature of the 
Ising and vortex unbinding transitions. 

It is worth noticing that an interesting realization of the FFXY model physics just 
outlined can be given in terms of a two dimensional JJA in a transverse magnetic field 
[4, 14, 27, 38]. Such a system is a periodic array of superconducting islands connected by 
Josephson links; each island is characterized through the modulus and the phase of the 
order parameter. The lowering of the temperature to its critical value forces each island 
to become superconductive but the whole array is in a resistive state as long as the phases 
have not acquired a long range order. The array reaches such a global coherence state at 
a lower temperature and its behaviour can be well described only in terms of the order 
parameters phases fi through the Hamiltonian: 

Hjj A = -Ej ^ cos (<Pi - <Pj) , (30) 

(ij) 

where Ej > is the Josephson coupling energy and the sum is over nearest neighbours. 
Let us notice the similarity of Eq. (30) with the Hamiltonian of a two dimensional XY 
model. In this picture it is known that the array becomes superconducting below the 
temperature Tj = ^rp-, where the phases are ordered and vortices may appear only in 
bounded pairs. Hence the ferromagnetic state for the XY model corresponds to the 
superconducting state for the JJA [13]. Such a regime is the classical one for a JJA, 
indeed the quantum fluctuations of the phases tpi are weak and the vortices are well 
defined objects, which form a Coulomb gas and behave as particles with masses. 

The analogy just outlined still works very well in the presence of a magnetic field 
perpendicular to the plane of the array, in this case a frustrated XY model can be realized; 
in particular, the full frustration f — \ corresponds to half a magnetic flux threading each 
cell. The ground state can be viewed as a pinned vortex lattice commensurate with the 
underlying periodic pinning potential leading to discrete symmetries, which add to the 
U (1) symmetry of the superconducting order parameter [4, 14]. Furthermore we must 
remark the striking role of fully frustrated JJAs as an experimental tool for the study of 
phase transitions in FFXY type models [33, 34]. 

Acknowledgments. We would like to thank E. Orignac for his careful reading and 
comments on a preliminary version of this work. G. N. was supported by a postdoctoral 
fellowship of the Ministere frangais delegue a l'Enseignement superieur et a la Recherche. 
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A TM on a torus 



The partition function of the TM model on the torus has the following factorized form 
(see [3]): 



Z(w c \t) = Z mr (w c \t)Zj(t), 



(31) 



in terms of the Moore- Read partition function Z MR (c = 3/2) and of the Ising partition 
function Zj (c = 1/2). They have the following expression: 

Z MR (w c \t) = \ X ^ c (0\w c \r)\ 2 + | X f c (0Mr)| 2 + |xf C (0K|r)| 2 , (32) 



^) = IXo(r)| 2 + U (V) 
The above characters are defined by: 



+ 



X-L{T) 



(33) 



Xo(r) 

x±S T ) 



1 V(r) + V 


l9 4 (0\r)\ 
V(r) J 


/»3(0|r) 


le 4 (o\r)\ 


/ V(r) V 





i e 2 (o\T) 

Mr) ' 



(34) 
(35) 
(36) 



which express the primary field content of the Ising model with Neveu-Schwartz (Z 2 
twisted) boundary conditions, and by: 



X ^ C (0\w c \t) = xo(t)K (w c \t)+ X i(t)K 2 (w c \t), 
xf C (0K|r) = X ±(r) (Mwclr) + K 3 (w c \t)) , 

16 

X 2 C (0\w c \r) = xi(r)K (w c \r)+xo(0\r)K 2 (w c \r) 



(37) 
(38) 
(39) 



which express the primary field content of the Z 2 -invariant c = | CFT. They are given 
in terms of a "charged" contribution (c = 1): 



K 2 i +i {w c \T) = —— e 

77 (r) 



2l+i 
4 



(2w c \At), V(/,*)e(0,l) 2 , 



(40) 



and a "isospin" one X/?( r ) ( c = 1/2), where w c = — : ln^ c is the torus variable of the 

2iri 

"charged" component. 
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